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Abstract 

We calculate penguin contributions to the lifetime splitting between the B s and the B d meson. In 
the Standard Model the penguin effects are found to be opposite in sign, but of similar magnitude 
as the contributions of the current-current operators, despite of the smallness of the penguin 
coefficients. We predict 

r(jBs) -l = (-1.2 ± 10.0) -10- 3 ' fBs > 



r(B d ) v ' ' V190MeV; ' 

where the error stems from hadronic uncertainties. Since penguin coefficients are sensitive to 
new physics and poorly tested experimentally, we analyze the possibility to extract them from 
a future precision measurement of r (B s ) jr [Bd) - Anticipating progress in the determination of 
the hadronic parameters e\, e 2 and fs a / fB d we find that the coefficient C4 can be extracted with 
an uncertainty of order | AC 4 | ~ 0.1 from the double ratio (t(B s ) — T(B d ))/ (t{B + ) — T(B d )), 
if lei — £2! is not too small. 
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1 INTRODUCTION 



1. Introduction 

The theoretical achievement of the Heavy Quark Expansion (HQE) [1] has helped a lot to un- 
derstand the inclusive properties of B-mesons. The measurements of lifetime differences among 
the b-flavoured hadrons test the HQE at the order (Aq C d / 'm b ) 3 '. Today's experimental infor- 
mation on the B-meson lifetimes is in agreement with the expectations from the HQE, but the 
present theoretical predictions still depend on 4 poorly known hadronic parameters B u B 2 , E\ 
and e 2 [2, 3] • Recently they have been obtained by QCD sum rules [4] . Lattice results are expected 
soon from the Rome group [5] and will allow for significantly improved theoretical predictions 
of the lifetime ratios. 

Weak decays are triggered by a hamiltonian of the form 



Here Qi and Q 2 are the familiar current-current operators, Q 3 . . . Q e are penguin operators and 
Qs is the chromomagnetic operator. Their precise definition is given below in (3). The factors 
Vckm and V' CKM represent the factors stemming from the Cabibbo-Kobayashi-Maskawa matrix 
and are specific to the flavour structure of the decay. Feynman diagrams in which the spectator 
quark participates in the weak decay amplitude induce differences among the various b-flavoured 
hadrons. Such non-spectator effects have been addressed first by Bigi et al. in [6] evaluating the 
matrix elements in the factorization approximation in which e x = e 2 = 0. Then Neubert and 
Sachrajda [2] have found that even small deviations of Ei,e 2 from zero drastically weaken the 
prediction of [6] for the lifetime ratio r(B + ) /r(B d ), which can sizeably differ from 1. On the 
other hand the deviation of r(B s ) / r(B d ) from unity has been estimated to be below 1 % in [2, 6] 
and the detailed analysis of Beneke, Buchalla and Dunietz [3]. Here r(B s ) is the average lifetime 
of the two CP-eigenstates of B s . 

Experimentally the ratio r(_B s )/r(_B d ) can also be addressed by the measurements of the corre- 
sponding semileptonic branching fractions. Since spectator effects in the semileptonic decay rate 
are negligible, one may use T(B s )/r(B d ) = B SL (B s )/B SL (B d ). 

So far only the effect of Qi and Q 2 has been considered in [2, 3, 6]. Taking into account the 
present experimental uncertainty and the fact that C\ and C 2 are much larger than C 3 _ 8 in the 
Standard Model this is justified. Yet once the lifetime ratio T(B s )/r(B d ) is measured to an accu- 
racy of a few permille, the situation will change: The smallness of \r(B s )/r(B l i) — 1 j is caused by 
the fact that the weak annihilation contribution of Q ij2 depicted in Fig. 1 almost yields the same 
contribution to the decay rates of B s and B d . The difference in the CKM-factors is negligible 
and the lifetime difference is induced by the small difference of the (c, c) vs. (c, u) phase space 
and by 577(3) p violations of the hadronic parameters. These effects suppress \t{B s ) / r(B d ) — 1 1 
by roughly an order of magnitude compared to \t(B + ) / r(B d ) — 1 1 . The contributions stemming 
from the penguin operators and the chromomagnetic operator, however, do not exhibit such a 
cancellation. Their contribution to the non-spectator rate of B s comes with the same power of 
the Wolfenstein parameter A = 0.22 as the contribution of Q\ y2 . In contrast the effects of Qz~s 
to the non-spectator rate of B d or B + are suppressed by two powers of A and are therefore neg- 
ligible. Hence one expects the contributions of Q3-6 and Q 8 to \r(B s )/T(B d ) — 1| to be of the 



H 




(1) 



Figure 1: Non-spectator (weak annihilation) 
contribution to the B s decay rate involving two 
current-current operators. The corresponding 
diagram for the B d decay is obtained by re- 
placing s by d and the upper c by u. 



Figure 2: Weak annihilation diagram involving 
one penguin operator Q^-q. Penguin contribu- 
tions to the non-spectator rate of the B d meson 
are CKM suppressed and therefore negligible. 



same order as those of Q 1 and Q 2 . T(B + )/r(B d ) is not modified, so that the phenomenologi- 
cal conclusions drawn from this ratio in [2] are unchanged. Observables sensitive to C 3 _ 8 like 
r(B s )/r(B d ) are phenomenologically highly welcome. The smallness of C 3 _ 8 is a special fea- 
ture of the helicity structure of the corresponding diagrams in the Standard Model. In many of its 
extensions the values of these coefficients can easily be much larger. Such an enhancement due 
to supersymmetric contributions has been discussed in [7]. Up to now the focus of the search for 
new physics has been on new contributions to C 8 [7]. Yet many interesting possible non-standard 
effects modify C 3 _ 6 rather than C 8 : New heavy particles mediating FCNC at tree-level or mod- 
ifications of the b-s-g chromoelectric formfactor affect C 3 _ 6 , but not C 8 . Likewise new heavy 
coloured particles yield extra contributions to C 3 _ 6 , e.g. in supersymmetry box diagrams with 
gluinos modify C 3 _ 6 . 

It is especially difficult to gain experimental information on the numerical values of the penguin 
coefficients C 3 _ 6 . Even penguin-induced decays to final states solely made of d and s quarks do 
not provide a clean environment to extract C 3 _ 6 : Any such decay also receives sizeable contri- 
butions from Q 2 via CKM-unsuppressed loop contributions [8,9]. In exclusive decay rates these 
"charming penguins" preclude the clean extraction of the effects of penguin operators [8]. In 
semi-inclusive decay rates like B — > X s $ the situation is expected to be similar. In inclusive 
decay rates such as the total charmless b decay rate the effect of "charming penguins" can be 
reliably calculated in perturbation theory. Yet these rates are much more sensitive to new physics 
contributions in C 8 rather than in C 3 _ 6 , because Q$ triggers the two-body decay b — > sg, while 
the effects of Q3-6 involve an integration over three-body phase space [9]. Notice from Fig. 2 
and Fig. 3, however, that this phase space suppression of the terms involving C 3 _6 is absent in 
the non- spectator diagrams inducing the lifetime differences. 

This work is organized as follows: In the following section we calculate the contributions to 
r(B s )/r(B d ) involving Q 3 _6 or Q 8 . Here we also obtain the dominant part of the radiative 
corrections to order a s . In sect. 3 we discuss the phenomenological consequences within the 
Standard Model and with respect to a potential enhancement of C 3 _ 8 by new physics. 



4 



2 PENGUIN CONTRIBUTIONS 



2. Penguin Contributions 



For the non-spectator contributions to the B s decay rate we need the \AB\ = \AS\ = 1- 
hamiltonian: 



with 



Qx = {sc)v„a ■ {cb)v^ A ■ 1, 

Q3 = (sb) v -A ■ {qq)v~A ■ 1, 

q=u,d,s,c,b 

Q5 = {sb) v ~A ■ {qq)v+A ■ 1, 



J'=l 



<?2 = (sc)y_A • (c6)y-A • 1 
q=u,d,s,c,b 

Qe = (s b )v-A ■ (qq)v+A ■ 5 

q=u,d,s,c,b 



(2) 



Qs = 



9 



8tt 2 



m b sa^(l + l5 )T a b-G%. 



(3) 



The colour singlet and non- singlet structure are indicated by 1 and 1 and V ± A is the Dirac 
structure. For more details see [9, 10]. In (2) we have set V u bV* s = 0(A 4 ) to zero. The diagram 
of Fig. 1 has been calculated in [2, 3] and yields contributions to the non-spectator part r non_spcc 
of the B s decay rate proportional to Cf , C\ ■ C 2 and C\. The result involves four hadronic matrix 
elements, which are parametrized by the B-factors B\, B 2 , e\ and e 2 [2]: 



•2 



fk M%B 2 



(B s 1^(1-75)667^(1-75)515 

(B s \s(l + l5 )bb(l- l5 )s\B. 

(^1^(1-75)^667^(1-75)^!^) = f B M\ 
(B s \s(l + l5 )T a bb(l- l5 )T a s\B, ' ! ' " ! 



1.^1 
-2 



= fk K 



B s e l 

r k^- 



(4) 



Here T a is the colour SU (3) generator, M Bs = 5369 ± 2 MeV and f Bs are the mass and decay 
constant of the B s meson. T(B s )/r(B d ) - 1 is proportional to r non - spcc (B d ) - T non - spcc (B s ). 
The main differences between the result of Fig. 1 for these two rates are due to the different 
mass of u and c and the difference between fs d and fs s - Hence the current-current parts of 
t(B s ) I r(B d ) — 1 proportional to Cf , C\ ■ C 2 or C\ are suppressed by a factor of z or A with 



mt 



0.085 ±0.023, 



1 - 



fl 



f 2 B ME 



0.23 ±0.11. 



(5) 



The result for A in (5) is the present world average of lattice calculations [11]. There are also 
SU(3)f violations in the B-factors, but they are expected to be small from the experience with 
those appearing in B° — i? -mixing. We want to achieve an accuracy of 2 permille in our pre- 
diction for T(B s )/r(B d ), which corresponds to an accuracy of 20-30% in T(B s )/r(B d ) — 1. 
Therefore we use the same Bi, B 2 , e\ and e 2 in t(B s ) and r(B d ). Likewise there is SU(3) F - 
breaking in the matrix elements of the b-quark kinetic energy operator and the chromomagnetic 
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Figure 3: Contribution of Q 8 to r non - spec (£ s ). 
In the Standard Model the diagram is of the 
same order of magnitude as radiative correc- 
tions to Fig. 2 and therefore negligible. Yet 
in models in which quark helicity flips occur 
in flavour-changing vertices | C 8 j can easily be 
ten times larger than in the Standard Model 
[7]. The contribution of Q 1 vanishes. 




Figure 4: Penguin diagram contribution to 
r non_spec (S s ). The final state corresponds to a 
cut through either of the (c, c) -loops. The con- 
tributions of Qi vanish by colour. This is the 
only NLO contribution to r (B s ) / T(B d ) — 1 in- 
volving Qi j2 without suppression factors of A 
or z. 



moment operator. These effects are suppressed by a factor of rrib/(AQCD • 167r 2 ) with respect to 
those discussed above. In [3] they have been estimated from heavy meson spectroscopy to be an 
effect of order one permille in t(B s )/t(B ( i) . 

We are now interested in the diagram of Fig. 2 involving one large coefficient C 1:2 and one small 
penguin coefficient C 3 _ 6 . Diagrams with two insertions of penguin operators yield smaller con- 
tributions proportional to Cf_ 6 and are neglected here. To order A 2 in H we have V' CKM = in 
(1) for the Bd system and penguin effects are only relevant in r(B s ). Hence the penguin contri- 
butions to T(B s )/T(B d ) — 1 do not suffer from the suppression factors z and A. Next we want 
to evaluate the diagram of Fig. 3 which encodes the interference of Qi,2 with the chromagnetic 
operator Q 8 . This part of r non_spec already belongs to the order a s and is small in the Standard 
Model, but it can be sizeable in the new physics scenarios discussed in [7]. 
We also must discuss radiative corrections to the contributions involving the large coefficients 
C\ and C 2 . Dressing the diagram in Fig. 1 with gluons gives contributions to r non ~ spcc for both 
B d and B s and therefore yield small corrections of order C 2 Aa s /n or less. The penguin diagram 
of Fig. 4, however, contributes only to Y non ~ spcc (B s ) m the order A 4 . Hence Fig. 4 yields an 
unsuppressed contribution of order C 2 a s /ir to t(B s ) /r(B d ) — 1 and cannot be neglected. The 
result of these penguin loop diagrams can easily be absorbed into the penguin coefficients 6*3-6 : 
In the result of the diagram of Fig. 2 one must simply replace Cj by 

C'- = Cf LO + ^C 2 Re [r 2j (l, V£, fi/m b )} , j = 3, . . . , 6. (6) 

Here r 2 j encodes the result of the penguin diagram and can be found in [9] in the NDR scheme. 
To cancel the scheme dependence of r we must also include the next-to-leading order (NLO) 
corrections to Cj as indicated in (6). More precisely: We must include the NLO mixing of C 2 
into Cj in C^ LO , j = 3, . . . , 6, but the penguin-penguin mixing only to the LO. The difference 
between these partial NLO-coefficients, which are tabulated in [9], and the full C^ LO, s has a 
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2 PENGUIN CONTRIBUTIONS 



j 


1 


2 


3 


4 


5 


6 


8 


Cf ] (fx = m b ) 
C'j(i_i = m b ) 


-0.249 


1.108 


0.011 
0.014 


-0.026 
-0.041 


0.008 
0.014 


-0.031 
-0.047 


-0.149 


Cf\^ = m b /2) 
C'jin = m b /2) 


-0.361 


1.169 


0.017 
0.017 


-0.036 
-0.045 


0.010 
0.015 


-0.048 
-0.058 


-0.166 


Cf\n = 2m b ) 
C' j ( f i = 2m b ) 


-0.167 


1.067 


0.007 
0.012 


-0.018 
-0.036 


0.005 
0.013 


-0.020 
-0.039 


-0.135 



Table 1: The effective Wilson coefficients C] denned in (6) for z = 0.085, a s (M z ) =0.118 and 
m b = 4.8 GeV. Varying z within the range given in (5) affects the Cj by 3-4 % and is negligible 

for our purposes, are the LO Wilson coefficients. 



negligible impact on our result. Here we bypass this technical aspect of scheme independence by 
tabulating the Cj's in Tab. 1. 

Our result for the non-spectator part of the B s decay rate reads: 

pnon-spec (£ g = _ ^| | |2 ^Izf B M Bs 8, + d 2 S 2 + b 1 B 1 + b 2 B 2 ] (7) 

1Z7T 



with 



ai 
a 2 

h 
b 2 



2C\ + AC 2 C'} [l-z] + \2zC 2 C' % + [1 + 2z\ —C 2 C 8 

J 7T 



= -[l + 2z] 



2C\ + AC 2 C' + — C 2 C 8 



71 



[C 2 + N c C l ]^l-z) 
- [1 + 2z\ [C 2 + N C C\ 



N, ^ Nr 









+ 6z 


) 







N, 



[C 2 + N c d] + 2 



a' 

3 N r 



(8) 



Here N c = 3 is the number of colours. By setting Cj, j = 3, . . . , 6, and C$ in (8) to zero one 
recovers the result of [2]. 3 The result for the non-spectator contributions to the B d decay rate 
reads [2]: 



■pnon— spec {^Q^ 



12tt 



with 4 



a( = 2Cl[l + -), 



\V cb V ud \ 2 (1 - zf f B M Bs (A - 1) \ai e 1 + a d 2 e 2 + bfB, + b d 2 B 2 ] (9) 



4 = -2C 2 2 (l + 2z), 



N, 



(c 2 + N c c l y(\ 



td 

u 2 



-j-(C 2 + N c C 1 ) 2 (l + 2z).(W) 

3 Notice that our notation of C\ and C 2 is opposite to the one in [2]. 

4 In the large N c limit one finds r non-spcc helicity suppressed in analogy to the leptonic decay rate. This shows 
that one cannot neglect the 0(1/N C ) terms. 
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When we combine (7-10) in order to predict r(B s )/r(B d ) — 1: 

t(B s ) i pnon-spec^^ _ pion-spec^ 



t(B, 



- 1 



K(z) 



ptotal 



2C\ ( £l - e 2 ) + 



+ O(10~ 3 ) 
(C. + A^d) 2 



— 3C 2 Z£i — — 



3 (C 2 + iV c d) 2 



iV, 



z5i 



(B 1 - B 2 ) 



C\ ( £l - 4s 2 ) + (C2+ o !! cC7l)2 (Si - 4B 2 ) 



2iV, 



+ 



4C 2 C; + (l + 2z) — C 2 C 8 



7T 



(ei - £2) 



+ 2 (C 2 + 7V c d) + (5! - S 2 ) 
- Az C 2 C' 4 (£1 + 2e 2 ) + 12* C 2 C£ £1 



(11a) 

(lib) 

(He) 
(lid) 
(He) 



+2z {C 2 + N C C X ) 
+ 0(2- 10~ 3 ) 



C' 



a 



-[C' 3 + jf)(B 1 + 2B 2 ) + 3(C' 5 + ^)B 1 



Here K(z) reads 
K(z) 



16n 2 \V ud \ 2 B SL 



mlh{z) [l + a s (fi)/(27r)h SL (V^)] 
0.060 



f B M Bs [1 - 2z] 



1-4(^-0.3) 



h. v 



0.105 V190MeVy 



(12) 
(13) 



In (12) we have used the common trick to evaluate the total width r total in terms of the semilep- 
tonic rate and the measured semileptonic branching ration B S l via r total = Y S l/B S l- fi and 
h S L are the phase space and QCD correction factor of T S l calculated in [12]. We use the no- 
tation of [9]. The approximation in (13) reproduces K(z) to an accuracy of 3%. The numerical 
value of hsh entering (13) corresponds to the use of the one-loop pole mass (~ 4.8 GeV) for m\,. 
For simplicity we have expanded K(z) and the terms in the curly braces in (11) up to the first 
order in z. The size of the error in (11) is estimated as 2 • 10 -3 . Its main source is the SU(3) F - 
breaking in the kinetic energy and chromomagnetic moment matrix elements appearing at order 
A 2 Q CD /ml of the HQE, which has been calculated to equal (0-1) • 10~ 3 in [3]. Then terms of or- 
der 16ir 2 AQ CD /mf can maximally be of the same order of magnitude. Conversely the remaining 
NLO correction of order C 2 A a s /-K and the CKM-suppressed contributions are much smaller. 
Likewise the £77(3) ^-breaking in e± : e 2 , B\ and B 2 is expected to be at the level of a few percent 
and therefore smaller than the present uncertainty in A. 

The first three lines (lla-llc) contain the result of the current-current operators calculated in 
[2, 3]. The remaining lines comprise the penguin effects. Note that the terms in (lld-lle) are 
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3 PHENOMENOLOGY 



neither suppressed by A nor by z. For z = the hadronic parameters in (1 1) only appear in the 
combinations e\ — e 2 and B x — B 2 , both of which are of order 1 /N c . The coefficients of B x — B 2 
suffer from numerical cancellations, e.g. 0.09 < C 2 + 3Ci < 0.57 (cf. Tab. 1), so that for most 
values of the input parameters only the terms involving e 1 and e 2 in (11a), (lib) and (lid) are 
important. 

Finally we discuss a potential systematic uncertainty: The derivation of (1 1) has assumed quark- 
hadron duality (QHD) for the sum over the final states. QHD means that inclusive observables 
are unaffected by the hadronization process of the quarks and gluons in the final state. The new 
results for inclusive observables in B decays presented at the 1997 summer conferences are 
consistent with QHD [13]. There are two potential sources of QHD violation in our problem: First 
it may be possible that the spectator decay rate of the 6-quark is affected by the hadronization 
process. Yet the ballpark of this effect is independent of the flavour of the spectator quark and 
cancels out in the ratio r(B s )/r(Bd). S77(3)ir-breaking can only appear in the hadronization of 
the final state antiquark which picks up the spectator quark and we do not expect the 577(3) p- 
breaking in the spectator decay rate to be larger than the SU (3) ^-breaking in the (AQ C£ )/m 6 ) 2 - 
terms of the HQE. This effect should further not depend on whether the hadron containing the 
spectator quark recoils against other hadrons or against a lepton pair. Hence one can control the 
S77(3) ^-breaking in the spectator decay rate by comparing the hadron energy in semileptonic B d 
and B s decays. More serious is a potential violation of QHD in the non-spectator contribution 
pnon-spec ^ se if j n a theoretical analysis for the similar case of the width difference Ar Bs of 
the two B s eigenstates the size of QHD violation has been estimated to be moderate, maximally 
of order 30%. We can incorporate this into (11) by assigning an additional error of ±0.3 to A. 
In any case the issue of QHD violation in lifetime differences will be experimentally tested in 
the forthcoming years, when high precision measurements of T(B + )/r(Bd) and of Ar^ s are 
confronted with accurate lattice results for the hadronic parameters. 

3. Phenomenology 

In the following we want to investigate the numerical importance of the penguin contribution. 
Then we analyze which accuracy is necessary to detect or constrain new physics contributions to 
C3-6 by a precision measurement of T(B s )/r(B d ). 

The three main hadronic parameters entering (11) are A, f Bs and E\ — e 2 , while Bi and B 2 
come with small coefficients. The canonical sizes of the B -factors are = 0(1/ N c ) and B^ = 
l + 0(l/N c ). An important constraint on the e/s is given by the measured value of t(B + ) / r(B d ) 
[2]. The result of [2] for r non - spcc (B+) is obtained from (9) by replacing the af, bfs with 





a u 2 = b u 2 = 0. (14) 



The experimental world average [14] 



t(B+) 



= 1.07 ±0.04 



(15) 



r(B d ) 
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A = 0.12 


A = 0.23 


A = 0.34 




-0.3 


-0.1 





-0.3 


-0.1 





-0.3 


-0.1 





e 2 = -0.3 


4.3 


1.9 


* 


4.3 


4.5 


* 


4.3 


7.1 


* 


e 2 = -0.1 


3.7 


1.3 


* 


1.1 


1.3 


* 


-1.5 


1.3 


* 


£2 = 0.1 


* 


0.6 


-0.6 


-2.2 


-2.0 


-1.9 


-7.4 


-4.6 


-3.2 


e 2 = 0.3 


* 


-0.1 


-1.3 


* 


-5.3 


-5.2 


-13.3 


-10.5 


-9.1 



Table 2: Standard Model prediction for 10 3 • [t(B, s )/t(B (1 ) - 1] obtained from (11) for f Bs = 
190 Me V, fi = m b = 4.8 GeV, z = 0.085, a s (M z ) = 0.118 and B x = B 2 = 1. The entries 
marked with * are in conflict with the experimental constraint (15), which also implies S\ < 0. 
There is an overall error of ±2.0 (see 1 1) for all entries. 



£i - £2 = 


-0.5 


-0.3 


-0.1 


0.1 


0.2 




peng 


cc 


peng 


cc 


peng 


cc 


peng 


cc 


peng 


cc 


f Bs = 160 MeV 


3.9 


-8.8 


2.3 


-4.9 


0.8 


-1.0 


-0.8 


2.8 


-1.5 


4.7 


f Bs = 190 MeV 


5.4 


-12.4 


3.3 


-6.9 


1.1 


-1.5 


-1.1 


4.0 


-2.2 


6.7 


f Bs = 220 MeV 


7.3 


-16.6 


4.4 


-9.3 


1.5 


-2.0 


-1.4 


5.3 


-2.9 


9.0 



Table 3: The columns labeled with 'peng' list the penguin contribution to 10 3 ■ [t(B s ) / r(B d ) — 1] 
as a function of E\ — e 2 and f Bs . The other input parameters have little impact on the size of the 
penguin contribution. The current-current part of 10 3 • [t(B 3 ) / r(B d ) — 1] is listed for e\ = —0.1 
and A = 0.23. For the remaining parameters see Tab. 2. 



leads to the following constraint: 

£1 , ( _ , ± o,)(^) 2 (^)W 2 + o.05. 

In [4] the £j's and Bi's have been calculated with QCD sum rules within the heavy quark effective 
theory (HQET). The results are e x {n = m b ) = -0.08 ± 0.02 and e 2 (fi = m b ) = -0.01 ± 0.03 
and Z?i )2 = 1 + 0(0.01). In view of the smallness of the £j's, however, it is conceivable that 
other neglected effects are numerically relevant.For example a NLO calculation of the matching 
between HQET and full QCD amplitudes replaces in (7) and (9) by £j + diB-i, where d { is a 
coefficient of order a s (m b )/ir. Here we will consider the range \ei\, \e 2 \ < 0.3, and further obey 
(15). 

In Tab. 2 we have tabulated t(B s ) / r{B d ) — 1 for various values of A and e\, e 2 . We have further 
split T(B s )/T(B d ) — 1 into its current-current part consisting of (1 la-1 lc) and the new penguin 
part involving C' 3 _ 6 , C 8 . These results can be found in Tab. 3. From Tab. 3 we realize that the 
penguin contributions calculated in this work are comparable in size, but opposite in sign to 
the current-current part obtained in [3]. This makes the experimental detection of any deviation 
of T{B s )/T(B d ) from 1 even more difficult, if the penguin coefficients are really dominated by 
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Standard Model physics. The results of Tab. 2 can be summarized as 




r(B s ) 



1 



(-1.2 ±8.0 ±2.0) • 10 




Here the first error stems from the uncertainty in e\ and e 2 and will be reduced once lattice 



results for the hadronic parameters are available. The second error summarizes the remaining 
uncertainties. If A and e 2 simultaneously aquire extreme values, r(B s )/r(B d ) — 1 can be slightly 
outside the range in (17) (see Tab. 2). 

Today we have little experimental information on the sizes of the penguin coefficients. Their 
smallness in the Standard Model allows for the possibility that they are dominated by new 
physics. The total charmless inclusive branching fraction Br(B — > no charm) is a candidate 
to detect new physics contributions to Cg [7], but it is much less sensitive to C 3 _6 [9]. The de- 
creasing experimental upper bounds on Br(B — > no charm) [14] therefore constrain C$ but 
leave room for a sizeable enhancement of C 3 _ 6 . Now (11) reveals that T(B s )/T(B d ) is a com- 
plementary observable mainly sensitive to C 4 , while C 8 is of minor importance. As mentioned 
in the introduction, many interesting new physics scenarios affect C 3 _ 6 , but not necessarily C 8 . 
We remark here that we constrain ourself to new physics scenarios, in which the CKM factors 
of the new contributions are the same as the ones of the Standard Model. This is fulfilled to a 
good approximation in most interesting models [7]. Now any new physics effect modifies C 3 _ 6 
at some high scale of the order of the new particle masses, while the Wilson coefficients entering 
(1 1) are evaluated at a low scale // m^. The renormalization group evolution down to /j, m b 
mixes the new contributions to C 3 _ 6 . New physics contributions AC 3 _ 6 (yU = 200GeV) affect 
C 4 (At = 4.8 GeV) by 

AC 4 (/x = 4.8 GeV) = -0.35 AC 3 (200 GeV) + 0.99 AC 4 (200 GeV) 

-0.03 AC 5 (200 GeV) - 0.22 AC 6 (200 GeV). 

Observe that AC 4 (200 GeV) = —0.05 already increases C' 4 (m,b) by more than a factor of two. 
Clearly the usefulness of r(B s )/T(B d ) to probe C 3 _ 6 crucially depends on the size of \e 1 — e 2 \ 
and f Bs . We now investigate the sensitivity of t(B s ) / r(B d ) to AC 4 (/x = m b ) in a possible future 
scenario for the hadronic parameters. We assume 

e x = -0.10 ±0.05, e 2 = 0.20 ±0.05, B U B 2 = 1.0 ±0.1, 

f Bs = (190 ±15) GeV, A = 0.23 ± 0.05, m b = (4.8 ± 0.1) GeV. (18) 

The assumed accuracy for f Bs will be achieved, once more experimental information on the B s 
system is obtained, e.g. after the detection of B s — £? s -mixing. Also a more precise measurement 
of fo 3 is helpful, because lattice QCD predicts the ratio /b s //d s much better than f Bs [11]. The 
error bars of the other hadronic parameters likewise appear within reach, if one keeps in mind that 
information on e\ and e 2 will not only be obtained from the lattice but also from other observables 
like T(B + )/r(B d ). Experimental progress in (15) and a next-to-leading order calculation of the 
coefficients in (14) and (10) will significantly improve the constraint in (16). In Fig. 5 we show 
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Figure 5: Dependence of r(S s )/r(S d ) — 1 on a new physics contribution AC 4 . The shaded area 
corresponds to a variation of the input parameters within the range of (18). The horizontal lines 
mark the Standard Model range corresponding to AC 4 = 0. 



the dependence of r(B s ) / r(B d ) — 1 on AC 4 (/z) for the scenario in (18). A cleaner observable is 
the double ratio 

r(B s ) - r(B d ) B SL (B S ) - B SL (B d ) 

r(B+)-r(B d ) B SL (B+) - B SL (B d y K 

which depends on Si, e 2 and A, while the dependence on f B and m b cancels. The corresponding 
plot for the parameter set of (18) can be found in Fig. 6 

We find a smaller error band for (t(B s ) - r(B d )) / (t(B + ) - r(B d )) than for T(B s )/r(B d ) - 1. 
If AC 4 < -0.075 or AC 4 > 0.140, we find the allowed range for (r(B s ) - t(bJ))/ (t(B + ) - 
r(B d )) incompatible with the Standard Model. An experimental lower bound r(B s )/r(B d ) > 
1.005 would indicate a new physics contribution AC 4 < —0.063 in our scenario. Likewise the 
experimental detection of a sizeable negative lifetime difference r(B s ) — r(B d ) may reveal non- 
standard contributions to C 4 of similar size as its Standard Model value. Fig. 6 shows that e.g. 
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Figure 6: Dependence of (t(B s ) — r(B d )) j \t(B + ) — r(B d )) on AC 4 for the parameter set in 
(18). This double ratio depends on f Bs and f Bd only through A, and the factor of m^ 3 in (11) 
cancels. 



the bound r(B s ) - r(B d ) < -0.20(r(S+) - r(B d )) would indicate AC 4 > 0.051. We conclude 
that the detection of new physics contributions to C 4 of order 0.1 is possible with precision 
measurements of r(B s )/r(B d ). 

4. Conclusions 

We have calculated the contributions of the penguin operators Q 3 -a, of the chromomagnetic 
operator Q 8 and of penguin diagrams with insertions of Q 2 to the lifetime splitting between the 
B s and B d meson. In the Standard Model the penguin effects are found to be roughly half as big 
as the contributions from the current-current operators Q 1 and Q 2 , despite of the smallness of the 
penguin coefficients. Yet they are opposite in sign, so that any deviation of t(B s ) — r(B d ) from 
zero is even harder to detect experimentally. Assuming a reasonable progress in the determination 
of the hadronic parameters a precision measurement of T(B s )/r(B d ) can be used to probe the 
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coefficient C 4 with an accuracy of |AC 4 | = 0.1. Hence new physics can only be detected, if 
C 4 is dominated by non-standard contributions. The sensitivity to C 4 depends crucially on the 
difference of the hadronic parameters e 1 and e 2 . For the extraction of C 4 the double ratio (r(B s ) — 
T(B d ))/(T(B + ) - r(B d )) turns out to be more useful than T(B s )/r(B d ). 
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